Pathways to faceting of vesicles 
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An interplay between geometry, topology, and order can lead to geometric frustration that pro- 
foundly affects the shape and structure of a curved surface. In this paper we discuss how such effects 
can result in faceting of elastic vesicles. We show that under the right conditions an assortment of 
regular and irregular polyhedral structures can be the low energy states of such elastic objects. In 
particular, we discuss how topological defects necessarily present in any crystalline lattice confined 
to spherical topology naturally lead to the formation of icosahedra in a homogeneous elastic vesicle. 
Furthermore, we show that introducing heterogeneities in the elastic properties or allowing for a 
non-linear bending response of a homogeneous system opens non-trivial pathways to the formation 
of faceted, yet non-icosahedral structures. 



INTRODUCTION 



While most of the lipid bilayer vesicles found 
in biological systems are smooth, this is not al- 
ways the case and one often observes faceted struc- 
tures. Prominent examples in the nature are pro- 
tein capsids of some viruses, bacterial organelles 
car boxy somes and square-shaped bacteria found 
in saturated brine pools. 2 On the experimental side 
faceted forms are found in small phosphatidylcholine 
vesicles below the freezing temperature^ some 
cat ionic vesicles^ cross-polymerized vesicles j£ vesi- 
cles assembled from a mixture of anionic and cationic 
amphiphiles, 5,7,8 vesicles formed from block copoly- 
mers with liquid-crystalline side chains j2r— giant 
fuller enesji 2 - and gold nanocages^ While very dif- 
ferent in their structure, all these systems share the 
common feature of possessing pathways to lower the 
total energy by forming stable facets. Over the past 
four decades a substantial research effort has been 
devoted to understanding the shapes of smooth vesi- 
cles. Only in the past couple of years there has 
been growing interest in developing theoretical un- 
derstanding of the faceted structures. 

To understand the origin of faceting it is essen- 
tial to understand the physics of smooth vesicles. 
Early in evolution it became clear that in order 
to construct complex life it was necessary to be 
able to hierarchically compartmentalize and sepa- 
rate different parts of an organism. The compart- 
mentalization was achieved with the aid of vesi- 
cles made of thin membranes. Biological mem- 
branes are typically lipid bilayer structures inter- 
penetrated with a variety of membrane and trans- 
membrane proteins^ They are usually very thin 
(~ 5nm) compared to their lateral dimensions (tens 
to hundreds of nanometers to microns) and thus 
quasi-two-dimensional. While being ubiquitous in 
living organisms, vesicular membranes are not lim- 
ited to biological systems. Nano- and micro-size con- 
tainers can be constructed in the laboratory using 



lipid molecules ^ 15 ! 16 block-copolymers^ 7 - and even in- 
organic compounds, as in the case of gold nanocages. 

While the detailed structure and function of bio- 
logical membranes is exceedingly complex and not 
fully understood, much insight into their physical 
properties can be gained by constructing simple 
mesoscopic models. In a seminal paper, Helfrichi^ 
argued that due to its quasi-two-dimensional nature 
the dominant deformation of a membrane is bend- 
ing. Therefore it can be qualitatively described in 
terms of changes of the membrane's curvature. He 
constructed a free energy functional that relates the 
bending energy to the integral of the square of the 
mean curvature, H, and to the Gaussian curvature, 
K, i.e., 
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where Hq is the spontaneous curvature and k and kq 
are material-dependent bending and Gaussian rigidi- 
ties. It is evident that the Helfrich free does not ac- 
count for any change in the energy due to stretching 
or shear deformations. In other words, it is suitable 
for describing fluid membranes where molecules are 
confined to the membrane surface but are free to lat- 
erally diffuse within it. Fluidity is important for bio- 
logical membranes as it ensures that small apertures 
formed by transmembrane transport or as a result of 
a damage can quickly be closed; it also allows struc- 
tural reorganizations within the membrane itself, as 
in the case, for example, of lipid rafts ^ 

Homogeneous fluid vesicles are typically smooth 
as molecules quickly rearrange themselves to remove 
any stress that would be generated by the formation 
of sharp corners and edges. Sharp structures rep- 
resent regions with diverging curvature 2 ^ and are 
therefore energetically costly. While in real sys- 
tems the actual divergence is cut off at some molec- 
ular length scale, the energy penalty for forming 
sharp turns is comparable to the molecular ener- 
gies and far exceeds any energy penalty for form- 
ing long- wavelength undulations that are~ ksT. In 
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order to form sharp structures, therefore, it is nec- 
essary to have a mechanism that would either re- 
duce the energy cost associated with the formation 
of bends or prevent the flow of the molecules or im- 
pose additional "in- vesicle" order (e.g., liquid crys- 
talline) that would dictate the underlying geome- 
try. The former can be achieved by asymmetry in 
the composition between two leaflets of the bilayer 
induced by a phase segregation of the two molec- 
ular species, thus resulting in local non-zero spon- 
taneous curvature 2 ^ or by segregating excess am- 
phiphiles along long ridges and inducing a spon- 
taneous curvature that is commensurate with the 
dihedral angle between two faces i2£ In this paper 
we discuss mechanisms that lead to faceting of solid 
(crystalline) vesicles. 

This paper is organized as follows. In Section Hi Al 
we discus pathways to the faceting of elastic (teth- 
ered) membranes that describe viral capsids or bi- 
layer systems cooled below the gelation transition or 
for which the molecules are cross-polymerized (i.e., 
tethered to each other). In particular, in Section 
III Al we show how topological defects can act as the 
seeds of buckling into icosahedra providing a sim- 
ple explanation for the observation that viruses with 
icosahedral symmetry take spherical shape if their 
radius is small but are well faceted icosahedra if they 
are larger in size. In Section |TTB] we show how the 
presence of elastic inhomogeneities also can lead to 
faceting into regular and irregular polyhedra differ- 
ent from the ubiquitous icosahedron. Finally, in Sec- 
tion [TTC] we argue that faceting can also occur in ho- 
mogeneous elastic vesicles if one allows for non-linear 
bending response in terms of a critical curvature. 



II. FACETING OF ELASTIC MEMBRANES 
A. Buckling into icosahedra 

One of the hallmarks of viruses is a regular struc- 
ture of the capsid. The unit building blocks of 
the capsid are capsomeres, robust protein com- 
plexes, « lOnm in diameter that tile the capsid in 
a crystalline array. Lattice structures of icosahe- 
dral viruses can be described in terms of pairs of 
non-negative integers, (p, q) that form a T— number, 
T = p 2 + q 2 + pq 23 The existence of a lattice com- 
bined with the spherical topology is at the center of 
the argument for size-driven buckling into icosahe- 
dra. In order to explain the mechanism behind this 
buckling phenomenon we first need to briefly review 
the complex question of the existence of long-range 
order in curved spaces. 

Constructing a regular lattice in a plane (or any 
other zero Gaussian curvature surface) is fairly easy. 



The simplest example is the ground state of identi- 
cal point-like particles in an external confining po- 
tential interacting via long range repulsive interac- 
tions, e.g., identical charges confined to a planar re- 
gion. The ground state is a triangular lattice with 
each particle having exactly six equidistant nearest 
neighbors and a lattice constant that is determined 
by the density. Such a state is not frustrated and 
is stress free. This is no longer the case when one 
tries to perform a similar construction on the surface 
of a sphere. Any attempt to map a regular lattice 
onto a surface of the sphere will result in an ex- 
cess of particles and their mutual overlap. In other 
words, in order to build smooth covering of a sphere 
it is necessary to remove parts of the original planar 
lattice. As a result, even in the ground state any 
crystalline lattice on a sphere will necessarily have a 
finite number of sites that have coordination differ- 
ent than six, i.e., so called topological defects. Start- 
ing from the Euler's polyhedral formula (after all a 
tiling of a sphere can be thought of as a polyhedron!) 
which states that for any polyhedron V — E + F = %, 
where V is the number of vertices, E is the number 
of edges, F is the number of faces, and x is the Euler 
characteristic of the polyhedron, one can show that 
the total number of vertices of a crystalline lattice 
on a sphere has to satisfy the following relation^ 
(6 — z) N z = 6% = 12. N z is the total number 
of sites with coordination z and we have explicitly 
used that for a sphere x = 2. The configuration with 
the lowest number of defects on a sphere is A^ 5 = 12, 
i.e., twelve five- fold coordinated sites embedded in a 
regular triangular lattice. We note that, somewhat 
counterintuitively, this configuration is not always 
the ground state of the system as shown by Bowick, 
et al.: 25 for reasonable defect core energies the to- 
tal energy can be lowered by forming pairs of five- 
and seven-fold defects that emanate from each of 
the twelve five- fold coordinated sites. In the follow- 
ing discussion we will ignore this complication and 
assume that only twelve five- fold defects are present. 

Despite their appearance five-fold disclination de- 
fects are not local since they owe their existence to 
removing whole sections of a planar lattice so that 
it conforms to the surface of a sphere. As such they 
are endowed with considerable elastic energy. To 
estimate how the elastic energy of a disclination de- 
pends on its size we assume that it is surrounded by 
a Hookean elastic medium. The elastic energy of a 
two-dimensional crystal is^ 6 - E e \ — f d 2 rA % i kl UijUku 
where (r) = - (diUj + djUi + dkUidkUj) is the 
strain tensor field with u\ (r) (i = 1,2) being com- 
ponents of the displacement field and A l ^ kl is the 
material-dependent elastic tensor. For a hexagonal 
lattice, the elastic tensor has only two independent 
components, 26 i.e., two Lame coefficients ji and A 
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and the elastic energy simplifies to2& 

Eei = \j d 2 r{\ (u u ) 2 + 2/*u?.) , (2) 

and we assume summation over pairs of repeated 
indices. 

A planar five-fold disclination defect can be con- 
structed by removing a 60° angular section of a tri- 
angular lattice and sewing the two cut edges to- 
gether. During such a deformation, all points move 
along an arc traversing a path of length ex r, where r 
is the radial distance from the disclination. In other 
words, the displacement field produced by a five-fold 
disclination has a magnitude \u(r)\ cx \f\ resulting 
in a constant strain field, Uij (r) = const. This ob- 
servation immediately leads to the estimate of the 
stretching energy, 

E^~R 2 , (3) 

where R is the disclination radius. A detailed 
calculation 27 finds E h ~ fold = (Y/32tt)R 2 , where Y 
is the two-dimensional Young's modulus. 

If the disclination is allowed to buckle out of its 
plane it can lower its energy by forming a conical 
structure with the apex at the disclination. In this 
case, the stretching energy is relieved at the expense 
of a bending penalty. The bending penalty is given 
tqp2Z>2& 

j dS [2kR 2 + k g K) . (4) 

If the boundary of the integration domain does 
not change during buckling, the Gauss-Bonnet 
theorem 2 ^ ensures that the Gaussian curvature term 
does not contribute to the energy and one needs to 
retain only the mean curvature contribution. For 
a cone one of the principal curvature directions is 
along its slant height (i.e., = 0) while the other 
is perpendicular to it (i.e., K2 cx 1/r, with r being 
the distance from the apex). Therefore, we immedi- 
ately obtain, H cx leading to the estimate for the 
bending energy 

*t 5 L~i°g(W, (5) 

where a is some microscopic length scale of the or- 
der of the lattice spacing. A detailed calculation^ 7 - 
finds 4 end ~ (7r /3)ttlog(fl/fl b ) + stYRl where 
R b « ^/IMk/Y. By comparing Eqs. §3§ and ([5]) 
it is clear that for small radii, the elastic stretch- 
ing energy will be smaller than bending energy and 
the disclination will remain flat at the expense of 
stretching its surrounding. As R increases, the log- 
arithm grows less rapidly than R 2 and it becomes 




(a) (b) 

Figure 1: Disclination-driven buckling into an icosahe- 
dron. The transition is triggered if the dimensionless 
FvK number 7 exceeds a critical value 7 ~ 10 2 . 30 Sim- 
ulated structures are constructed using (p,q) = (12,3) 
T— number triangulation with « 2 x 10 3 vertices and ra- 
dius R w 11.5/o, with lo = 1 being the average lattice 
spacing, (a) 7 « 170 and (b) 7 « 3500. Red dots desig- 
nate five- fold defects. All snapshots are generated with 
the Visual Molecular Dynamics (VMD) package 3 -^ and 
rendered with the Tachyon ray-tracer.— 



energetically favorable to buckle into a cone. A de- 
tailed calculation 27 shows that the buckling transi- 
tion occurs if the dimensionless ratio, the so called 
Foppl-von Karman or FvK number, 7 = > 154. 

Another consequence of their non-local character 
is that disclination defects, much like charges, ex- 
ert long-range forces on each other. This long-range 
force is transmitted via elastic deformation of the 
medium with the speed of sound playing the role 
of the speed of light in electrostatics. Disclinations 
with the same topological charge (s = 1 for a five- 
fold defect) repel each other. Therefore, on a sphere 
disclinations take a configuration that maximizes 
their separation, which is an icosahedron for twelve 
defects^ 

Moving from the planar to spherical geometry 
does not change the functional forms of the stretch- 
ing and bending energies, albeit renormalizes nu- 
merical prefactors due to partial screening of the 
strain induced by the non-zero Gaussian curvature 
of the sphere. 25 As in the planar case, therefore, 
the onset of buckling solely depends on the relative 
strengths of stretching and bending energies. Since 
the stretching energy rapidly grows with the radius 
of the sphere once a critical value is reached all de- 
fects simultaneously buckle into cones thus forming 
an icosahedral vesicle. 30 In Fig. [1] we show the shape 
of the shell below and above the buckling transition. 
The transition is not sharp but a rounded smooth 
crossover, with vesicles having the FvK number 7 < 
150 being spherical, for 200 < 7 < 1500 being no- 
ticeably buckled, and for 7 > 2000 being very sharp. 
For extremely large FvK numbers (7 > 10 7 ) the vesi- 
cle is non-extensible and one observes a characteris- 
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tic ridge-like structures with very interesting scaling 
properties, e.g., E ridge /tz ~ a 7/s (YL 2 /k) 1/6 , where 
a is the angle of the ridge and L is its length i 33 i 34 

Finally, we note that imposing internal pressure 3 - 5 - 
or a volume constrain t 35 ^ 6 does not qualitatively af- 
fect the buckling transition but shifts its onset to 
higher values of 7. 




B. Buckling of multicomponent vesicles 

An implicit assumption in the discussion thus 
far has been that the vesicle is made of an elas- 
tically homogenous material. Introducing hetero- 
geneities by allowing for spatially varying Young's 
modulus and bending rigidity can lead to an al- 
ternative buckling mechanism. For example, pres- 
ence of heterogeneities is believed to be responsi- 
ble for the structure of recently synthesized assem- 
blies of oppositely charged amphiphiles. 6,8 Observed 
structures are ~ lOOnm size bilayer vesicles that are 
faceted but not icosahedral. SAXS/WAXS measure- 
ments combined with atomistic and semi-atomistic 
molecular dynamics simulations have revealed the 
existence of strongly correlated crystalline domains 
within the bilayer. 8 A typical domain size was found 
to be ~ 25nm. Simulations show that the crystalline 
orders of two neighboring domains are not mutu- 
ally aligned and that the domains are separated by 
narrow regions of disordered, or even liquid, am- 
phiphiles. It is reasonable to expect that the bi- 
layer will be softer along those disordered boundary 
regions. These regions represent grain boundary de- 
fects with very different physical properties from the 
topological defects discussed in the previous section. 

In order to model a vesicle with inhomogeneous 
elastic parameters we generalize the expressions for 
stretching (Eq. (|2J)) and bending (Eq. (jlj) energies 
and allow the two Lame coefficients and the bending 
rigidity to depend on the spatial position, i.e., A = 
A(r), \i = /i(r), and n = k{t). For simplicity we 
assume that each of these coefficients can take only 
two values, denoted as "hard" and "soft" and loosely 
corresponding to strongly correlated facets and grain 
boundary regions, respectively. 

The elasticity theory of thin objects is non-linear 37 
and if one allows for spatially varying elastic param- 
eters (even if they have a simple binary distribution) 
it quickly becomes analytically intractable. It is nec- 
essary to devise a scheme that allows us to obtain 
low energy states of the system using numerical sim- 
ulations. The vesicle is represented as a discrete tri- 
angulated surface. In the case of modeling viral cap- 
sids, the vertices of the triangulation correspond to 
individual capsomeres. If the constitutive units are 
smaller, as is the case for amphiphilic vesicles, the 
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(c) (d) 

Figure 2: Snapshots of faceted structures with 3% (a), 
4% (b), 4.5% (c), and 5% (d) of the soft component (yel- 
low). Twelve five- fold dislocation defects are present as 
required by the spherical topology but appear to have no 
influence on the position of the edges and size of facets. 

triangulation can be thought of as a coarse-graining 
scheme that assigns discrete elements to the patches 
of actual molecules. Each patch is large enough to 
contain a sufficient number of microscopic degrees of 
freedom that molecular details are of no importance 
but small enough to be considered homogeneous. 

On a triangular lattice the simplest form of a dis- 
crete stretching energy for a multicomponent sys- 
tem can be obtained if we assume that the neigh- 
boring sites are connected to each other with har- 
monic springs with spatially-dependent spring con- 
stants 6ij and a uniform rest length Iq 27 

(iJ) 

where Uj = \fi — fj\ is the Euclidean distance be- 
tween sites at fi and fj : Sij G {ehard, £soft}, and the 
sum is carried out over all pairs of nearest neighbors. 
It is convenient to set Iq = 1 and use it as a unit of 
length. 

A proper discretization of the bending energy 
is more complicated with a number of approaches 
proposed in physics, applied mathematics, and 
computer graphics literature £Z££r— Different ap- 
proaches vary in the level of computational complex- 
ity and accuracy. For the present discussion it is suf- 
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ficient to adopt the simple discretization suggested 
by Seung and Nelson^ The discrete bending energy 
is computed as 

Ebend = \ ^2 k u (Hi - njf 

= ^S /<7 (1- cos (?) 

(i, J) 

where 0jj is the angle between unit-length normals 
nj and nj to the neighboring triangles I and J, and 
kjj is the position dependent discrete bending rigid- 
ity also picked from a binary set {nhard-, Ksoft}- The 
spring constant and the discrete bending rigidity are 
related to the Young's modulus and bending rigid- 
ity of the continuum theory by Y = 2e/V2 27 (with 
Poisson's ratio v = 1/3) and k = V3k/2 27 ^ 2 

Unlike single component systems where only 
shape changes are allowed, thus rendering standard 
numerical minimization techniques (e.g., conjugate- 
gradient minimization) applicable, for multicompo- 
nent vesicles it would be beneficial to have a strategy 
that simultaneously optimizes shape and the com- 
ponent arrangement. This can be achieved by using 
stochastic minimization techniques, e.g., the simu- 
lated annealing Monte Carlo method. Two compo- 
nents are assigned at random to vertices (or equiv- 
alently to edges or triangles) of the initial sphere 
such that the total fraction of the soft component is 
/. A Monte Carlo move consists of two stages: 1) a 
vertex is displaced by a vector Ar with components 
chosen at random from a uniform distribution in an 
interval [— C^o,C^o] (typically, £ « 0.05) followed by 
2) swap of types (soft^hard) of a pair of randomly 
selected vertices. In both stages moves are accepted 
according the Metropolis rules. During a simula- 
tion, various cooling protocols can be applied, e.g. 
linear, exponential, power-law, etc. We note that 
the annealing temperature is not the actual physi- 
cal temperature (characteristic energy scales in real 
vesicles significantly exceed fc#T and thermal fluc- 
tuations can be ignored; however an interesting ex- 
ample where this is not the case has recently been 
discussed 43 ) but rather a parameter that controls 
the rate at which Monte Carlo moves that increase 
energy are accepted. During the component swap 
stage the total fraction of each component is kept 
fixed, i.e., no moves that would transform type of a 
vertex, e.g., from being hard to being soft, are al- 
lowed. Finally, the component swap move does not 
correspond to an actual rearrangement of the mate- 
rial within the vesicle - the vesicle is assumed to be 
solid and no diffusion is allowed, but is a convenient 
simulation tool that allows sampling of the space 
of component arrangements. Throughout a simula- 
tion the triangulation is preserved, i.e., no edge-flip 



moves 44 that would create or annihilate defects are 
allowed. 

The main drawback of the simulated annealing ap- 
proach is that it typically converges very slowly once 
a minimum has been approached (for studied system 
sizes a minimum is typically reached after 10 5 — 10 6 
Monte Carlo sweeps) and there is no guarantee that 
the obtained structure is true a global minimum, 
i.e., the ground state configuration. Therefore, we 
refer to the obtained shapes as typical. Ultimately, 
there is no reason to expect that the experimentally 
observed shapes correspond to the equilibrium struc- 
tures. The realized shape will depend on the details 
of the assembly process in ways that are at present 
not understood. 

A detailed study^ 5 - of the shapes of two component 
elastic vesicles has found a wide variety of regular 
and irregular polyhedral shapes, some of which are 
shown in Fig. [2j If a small amount of soft component 
is added to an otherwise hard vesicle it arranges into 
branching lines. As the amount of the soft compo- 
nent is increased, these lines start to merge and form 
facets. Facets are nearly flat with very sharp bends 
along the soft ridges. The total number of facets in- 
creases with the amount of the soft component un- 
til a non-universal, size-dependent concentration is 
reached. Beyond this size-dependent concentration 
of the soft component it is no longer favorable to 
generate new facets but instead the soft component 
is randomly distributed inside the existing ones. Al- 
though each configuration still contains twelve five- 
fold defects, it appears that the size of the facets 
and the position of the soft ridges is not sensitive to 
their presence (Fig.[2])- We point out that it is possi- 
ble to entirely remove the residual strain induced by 
the defects either by appropriately tuning the spring 
rest lengths 46 or by directly working with reference 
and realized metrics j 47 i 48 Such an approach 49 shows 
that indeed faceting of the soft-hard two component 
vesicles occurs even if the defect contribution to the 
elastic energy is removed. 

Thus far we have not assumed any mixing penalty 
for the two components. The components segre- 
gated solely as a result of the bending preference. 
If a mixing penalty, i.e., an effective line tension, T, 
is introduced it competes with the organization of 
the components that is preferred by elasticity and 
leads to an even richer variety of shapes. 50 For a 
weak line tension elasticity dominates and the soft 
component forms highly bent ridges separating flat 
hard faces. The width of these soft ridges grows as 
the line tension increases. If V is increased even fur- 
ther, the faces start to coarsen and merge together, 
until the components fully separate at Tlo/hZhard ~ 
2Q-2 50 g U gg es ting that even a small amount of mix- 
ing penalty can lead to full phase segregation. A 
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(c) (d) 

Figure 3: Snapshots of faced elastic vesicles with various 
mixing penalties between soft (yellow) and hard (gray) 
components for a system with an equal amount of hard 
and soft vertices. Mixing penalty increases from (a) to 

(d). 




Figure 4: Bending energy (in units of the bending rigid- 
ity k) as a function of the angle between unit-length 
normals to two neighboring triangles without a critical 
angle (black), with C — ^ (red), with C = ^ (green), 
and with 9 C = f (blue). 



subset of shapes of two-component vesicles with line 
tension is shown in Fig. [3j 

We note that line tension can lead to very inter- 
esting shapes in two— and three-component liquid 
vesicles. 52 



C. Faceting due to critical curvature 

Some of the faceted shapes observed in experi- 
ments as well as those obtained in simulations (e.g., 



Fig. [2]) show sharp ridges with a large dihedral an- 
gle between neighboring facets. Once the radius of 
curvature becomes comparable with the molecular 
length scales the deformation is no longer locally 
small and microscopic details of the vesicle start to 
play a role. Eq. © was derived 27 under the as- 
sumption that curvature is smooth and slowly vary- 
ing between neighboring points, such the bending 
energy could be expanded in terms of small devia- 
tions of the normal vectors to a planar surface and 
only leading terms retained. For a locally large de- 
formation this is no longer justifiable and the expres- 
sion for the bending energy needs to be modified. 
This is clearly a very hard problem that requires a 
detailed knowledge of the microscopic structure of 
the vesicle wall. Instead of constructing such a de- 
tailed model, we make a simple assumption that be- 
yond a critical angle bending energy saturates. In 
other words, Eq. ([7j) is modified to 

ES^i = R^2(l-F(0 I j)), (8) 

where 

F{e) = [ cob (J) for Z0 e ( 

v 7 [ cos (0 C ) for > C v 7 

and as in Eq. ([7]) 9jj measures the angle between 
normals to a pair of neighboring triangles and C is 
the critical value of that angle beyond which the 
bending penalty saturates. Plot of the E^^d vs - 
is shown in Fig. 2J In this model C is treated 
as an input parameter. In essence, 6 C controls the 
deformation angle beyond which the vesicle wall is 
sufficiently deformed that atomistic details start to 
be important. We assume that the vesicle is ho- 
mogenous sphere with uniform values for R and e, 
and a standard triangular lattice with twelve five- 
fold defects positioned at the corners of an inscribed 
icosahedron. While one could easily argue that the 
stretching energy would also need to be modified, 
for simplicity, we use the expression in Eq. (j6j) 
and note that test simulations show that the results 
are not qualitatively affected by this simplification. 
As before, the low energy states were found using 
simulated annealing Monte Carlo simulations. No 
constraints on volume or area were imposed. Mini- 
mum energy configurations were found within « 10 6 
Monte Carlo sweeps using a linear cooling protocol. 

We studied vesicles with « 6- 10 3 vertices ((p, q) = 
(18, 10)- chiral; (25,0)- achiral) over a range of val- 
ues of the FvK number (7) and critical angle (0 C ). 
In Fig. [5] we show a phase diagram of this system 
mapped with « 200 independent pairs of 7 and C . 
Five distinct regimes are identified. For low values of 
7 and moderate to large C (green region, denoted as 
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K) 2 y 

Figure 5: Phase diagram for a vesicle with critical bend- 
ing angle as a function of the FvK number (7) and the 
value of the critical angle (0 C , shown in degrees). We 
identify five distinct regimes. In the "sphere" (green) 
and "smooth icosahedra" (yellow) regimes C is large and 
does not affect the shape. For small 7 the vesicle is spher- 
ical and crosses over into a smooth icosahedron as 7 is 
increased. At lower values of C the presence of critical 
angle starts to play a role and the vesicle shape changes 
substantially. For small 7 elastic energy is completely 
dominated by the bending contribution and the vesicle 
takes shapes of irregular polyhedra (red region). For 
larger values of 7 stretching energy caused by topologi- 
cal defects cannot be ignored and the icosahedral sym- 
metry reemerges, but accompanied by sharp ridges (blue 
region). Finally, for very small values of 7 the vesicle is 
unstable and crumples (gray). 

"sphere" in Fig. [SJ bending penalty is too large and 
the vesicle prefers to remain spherical. In this re- 
gions, angles between normals on any pair of neigh- 
boring triangles are smaller than C and the exis- 
tence of an angle cutoff plays no role on the shape. 
In other words, we recover the situation discussed in 
Section Hi Al If 7 is increased while keeping C large, 
one crosses over into the regime where icosahedral 
shapes are favorable. Since no deformation exceeds 
C the transition occurs at the same value of 7 as in 
the case without the critical bending angle (yellow 
region, denoted as "smooth icosahedra" in Fig. [5]). If 
6 C is reduced the effects of the angle cutoff start to 
emerge affecting the low energy shapes. In the re- 
gion where 7 is small (red region denoted as "polyhe- 
dra" in Fig. [5]) one observes irregular polyhedra with 
sizes and positions of facets and edges that are insen- 
sitive to the presence of the five- fold defects. This is 
to be expected as for small 7 bending penalty com- 



pletely dominates over stretching and defect-driven 
buckling into a cone is suppressed. In this region 
faceting mechanism is similar to the buckling as the 
soft component discussed in the previous section. As 
7 is increased the topological defects start to play the 
role and the system transitions in the regime where 
they drive buckling. With a relatively small value of 
C , however, the system can further lower its energy 
by forming sharp edges that emanate from the de- 
fects. This region is denoted as "sharp icosahedra" 
(blue in Fig. [5]). We note that the direction of the 
sharp edges appears to depend on the details of the 
triangulation. Finally, if C is too small the vesicle 
becomes unstable and crumples (gray region denoted 
as "crumpled" in Fig. [5j). In this region, the simu- 
lations often fail to converge and instead get locked 
in high-energy states that often occur in the early 
stages of a simulation. In Fig. [6] we show typical low 
energy configurations in each of the five regimes. 



III. SUMMARY AND CONCLUSIONS 

We have showed that faceted structures can be 
low energy configurations of elastic vesicles as a 
result of a number of different mechanisms. In- 
compatibility between crystalline order and spher- 
ical topology leads to the presence of topological 
defects resulting in a residual stress even in the 
ground state. The stress can be in part relieved 
by buckling into cones seeded at the defects. Due 
to the high symmetry of the ground state distri- 
bution of the defects, the buckled structure is an 
icosahedron. This mechanism can qualitatively ac- 
count for the observation that smaller icosahedral 
viruses (e.g. Polyomaviruses—) are smooth (spheri- 
cal) while larger viruses (e.g. Mimivirus^) are buck- 
led. For viruses the size of each capsomer is suffi- 
ciently large (~ lOnm) that the capsid is a nearly 
perfect lattice with ~ 10 2 — 10 3 units and no addi- 
tional defects beyond the minimum twelve five-fold 
sites required by the topology. In case of lipid bilayer 
vesicles lipid molecules are much smaller (occupying 
an area ;$ O.bnm 2 ) and each vesicle is assembled 
from hundreds of thousands of lipids. If a crystal 
is formed on such a vesicle as a result of cooling or 
strong electrostatic correlations it is much harder to 
form a perfect lattice and one instead expects grain 
boundaries to form. This effect is even enhanced by 
the spherical topology that hinders formation of a 
uniform crystal to begin with. In the grain bound- 
ary region molecules are disordered and may even 
be in the liquid state. As a result, the wall is ex- 
pected to be softer. We have shown that if one 
allows for the presence of softer component in an 
otherwise crystalline vesicle, the low energy config- 
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(b) 



uration can be regular and irregular polyhedra that 
are not icosahedra. The soft component forms ridges 
that are highly bent and separate hard fiat facets. 
The faceting mechanism is distinctly different from 
the buckling into an icosahedron seeded at the five- 
fold sites. Furthermore, it is reasonable to expect 
that the bending response of the vesicle wall is not 
linear and significantly changes if the local deforma- 
tion becomes large. We introduced a simple model 
that takes into account such effects at the most ba- 
sic level. Even with a such simplified model we were 
able to show that a non-linearity in the bending re- 
sponse can drastically affect the low energy shapes 
and lead to interesting faceted structures. 



(d) 



(e) 

Figure 6: Snapshots of the low energy configurations of 
a vesicle constructed using (p, q) — (18,10) T— number 
triangulation with w 6 x 10 3 vertices and radius R ~ 
20.5/o • (a) For low a FvK number 7 « 50 and a large 
critical angle C — 20° the vesicle remains spherical, (b) 
For 7 « 10 3 and C 20° one observes a smooth buckled 
structure with icosahedral symmetry akin to the shapes 
discussed in Section Hi Al (c) For small 7 w 60 and small 
^ c = 5° we find an irregular polyhedral vesicle with facets 
being insensitive to the position of the defects, (d) As 
FvK number if increased to 7 ~ 10 3 and the critical 
angle set to C = 14° the structure is icoshahedral but 
with sharp ridges emanating for each defect, (e) Finally, 
at 7 ~ 2.5 • 10 3 and very small critical angle, C = 2° 
the vesicle is no longer stable and crumples. Red points 
indicate positions of twelve five-fold defects. 



Understanding the mechanisms that lead to 
faceting of vesicles goes beyond the realm of an aca- 
demic exercise and could have a far-reaching impact 
on nano-technology. For example, flat faces could be 
easier to functionalize or otherwise chemically treat. 
Biochemical reactions are expected to have different 
rates near flat compared to curved surfaces, which 
could be utilized to construct highly sensitive trans- 
port agents. With the ability to accurately control 
the shape of vesicles it could be possible to engineer 
vesicles that would, e.g., selectively target only par- 
ticular agents. We hope this work will broaden the 
interest in this emergent subject of modern nano- 
science. 
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